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A  theoretical  analysis  is  given  for  one-dimensional  laser  cooling  below  the  Doppler  limit  of  J  =  4 
ground-state  atoms.  The  laser  field  consists  of  a  pair  of  connterpropagating,  linearly  polarized,  low- 
power  beams,  whose  polarization  directions  differ  by  an  angle  6  (0  <  0  <  jr/2).  For  8  <C  1,  the 
effective  optical-pumping  time  is  shown  to  increase  strongly  near  the  nodes  of  the  standing  wave, 
and  the  cooling  force  can  be  much  larger  than  that  for  9  —  t/2.  Moreover,  for  6  <  1,  it  can  be 
shown  that  the  stimulated  p  rt  of  the  atomic  diffusion  is  reduced  considerably  as  compared  with 
that  for  9  =  jr/2.  As  a  consequence  it  is  possible  to  achieve  an  equilibrium  atomic  distribution 
that,  for  9  •<  1,  is  characterized  by  a  mean  kinetic  energy  that  b  lower  than  that  predicted  to 
occur  for  9  =  ir/2.  The  equilibrium  velocity  distribution  is  not  necessarily  Maxwellian,  and  thus  the 
temperature  of  the  atomic  ensemble  may  not  be  well  defined.  The  achievable  kinetic  energy  is  so 
small  that  the  cooled  atoms  may  be  trapped  in  the  vicinity  of  the  laser-field  nodes.  _  i^C  l 
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I.  INTRODUCTION 

I11  the  past  several  years  it  has  been  demonstrated  both 
experimentally  and  theoretically  that  the  laser  cooling  of 
neutral  atoms  can  lead  to  an  equilibrium  atomic  kinetic 
energy  Ek  much  smaller  than  that  given  by  the  Doppler 
limit  for  two-level  atoms,  that  is  Et  C  /iT,  where  P  is 
the  spontaneous  decay  rate  of  an  excited  atomic  state.  It 
has  been  shown  that  the  sub-Dopplcr-limit  laser  cooling 
is  directly  related  to  effects  arising  from  the  multilevel 
structure  of  the  atomic  ground  state  (Refs.  [1-7]).  Effi¬ 
cient  sub-Doppler-limit  cooling  is  obtained  using  a  very 
weak  laser  field  satisfying 

« f>-  +  r'-\  (i.i) 

where  /  =  E/t/h  is  the  llabi  frequency,  and  /1  the  dipole 
moment  of  the  atomic  transition.  The  laser  field  has  an 
amplitude  E  and  a  frequency  w/,  detuned  from  the  atomic 
resonant  frequency  us, \  by  an  amount  8  given  by 

6=usr,-usA.  (1.2) 

Although  experiments  arc  carried  out  for  three-  [8-1 1] 
and  two-dimensional  [12]  field  configurations,  detailed 
analytical  calculations  that  have  been  performed  for  one- 
dimensional  cooling  [1,2,4]  (an  .atom  is  driven  by  a  pair 
of  two  connterpropagating  laser  beams)  give  a  fairly  ac¬ 
curate  description  of  the  major  features  of  the  cooling 
process.  It  has  been  shown  that  sub-Doppler-limit  laser 
cooling  is  produced  in  the  so-called  linXIin  configuration, 
when  the  laser  field  consists  of  a  pair  of  countorpropa- 
gal.ing  beams  having  orthogonal  polarizations.  For  this 
configuration,  the  polarization  gradient  of  the  laser  field 
plays  a  central  role  in  the  cooling  process.  'Phis  gra¬ 
dient  results  in  different  optical  pumping  rates  and  ac 
Stark  shifts  of  the  ground-state  sublevels,  both  of  which 
arc  spatially  modulated.  When  an  atom  moves  with  a 
nonzero  velocity  along  the  laser-beam  direct  ion,  it  can 


be  optically  pumped,  on  average,  from  ac-Stark-shifted 
levels  of  lower  energy  to  ac-Stark-shifted  levels  of  higher 
energy,  provided  the  sign  jf  the  atom-field  detuning  is 
chosen  properly  (so-called  Sisyphus  efTect  [1]).  This  op¬ 
tical  pumping  results  in  a  damping  of  the  atomic  veloc¬ 
ity.  In  general,  the  larger  the  spatial  modulation  of  the 
difference  in  ac  Stark  shifts  between  the  levels  and  the 
longer  the  optical  pumping  time,  the  stronger  the  fric¬ 
tional  damping  force.  For  weak  fields,  the  optical  pump¬ 
ing  time  Tp  can  be  made  much  larger  than  the  excited- 
.  state  lifetime  T-1  resulting  in  a  friction  force  that,  for 
slow  atoms,  is  much  stronger  than  that  in  the  case  of 
Doppler  cooling. 

In  the  case  of  parallel  linear  polarizations  for  the  fields, 
there  is  no  field  polarization  gradient  and  the  friction 
force  is  similar  to  that  ror  Doppler  cooling.  One  might 
think  that  as  the  rotati  n  angle  0  between  the  field’s  po¬ 
larization  vectors  varies  from  0  to  zr/2,  the  effectiveness 
of  sub-Doppler-limit  cooling  would  gradually  increase, 
achieving  its  maximum  for  0  =  ir/2.  This  conclusion  is 
not  necessarily  correct.  For  small  angles,  (1<  1,  the  in¬ 
tensity  of  the  laser  field  varies  considerably  in  space,  and 
the  optical  pumping  time  strongly  depends  on  the  p 
tion  of  the  atom.  A  though  the  ac  Stark  shifts  diffe.  only 
slightly  for  small  0,  the  fact  that  the  optical-pumping 
time  is  increased  significantly  near  regions  of  low  field  in¬ 
tensity  can  lead  to  a  cooling  force  much  stronger  than 
that  corresponding  to  0  —  it/2.  Moreover,  owing  to 
the  decrease  in  the  difference  of  ac  Stark  shifts,  the  mo¬ 
mentum  diffusion  coefficient  that  characterizes  stochastic 
heating  of  atoms  is  much  smaller  for  small  angles  than  for 
0  =  7t/2.  The  increase  of  the  force  and  decrease  of  the  dif¬ 
fusion  may  eventually  lead  to  a  lower  equilibrium  atomic 
kinetic  energy  for  smaller  angles,  than  for  0  =  for  a 
range  of  field  intensities.  It  should  b"  noted,  however, 
that,  if  the  field  intensities  are  chose1  to  minimize  the 
atomic  kinetic  energy,  optimnl  cooling  >ccurs  for  0  =  4 
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It  is  the  purpose  of  this  paper  to  derive  equations  de¬ 
scribing  laser  cooling  of  atoms  having  a  J  =  \  ground 
state  for  arbitrary  relative  polarization  angle  0.  In  de¬ 
riving  the  results,  one  finds  that  standard  computational 
techniques  involving  continued  fractions  or  iterative  so¬ 
lutions  are  not  practical  for  a  range  of  atomic  velocities. 
An  alternative  analytical  approach  is  used  to  solve  for 
the  atomic-state  density-matrix  elements. 

In  Sec.  II  a  qualitative  picture  of  sub-Doppler-limit 
cooling  is  given  for  an  arbitrary  relative  polarization  an¬ 
gle  0  between  the  fields.  The  Fokker-Planck  equation  for 
the  density  matrix  in  the  Wigner  representation  is  de¬ 
rived  in  Sec.  III.  The  results  for  the  force,  diffusion,  and 
achievable  temperatures  are  discussed  in  Sec.  IV.  Conse¬ 
quences  and  implications  of  the  results  obtained  in  this 
paper  and  related  phenomena  are  discussed  in  Sec.  V. 

II.  QUALITATIVE  PICTURE 

The  laser  field  is  represented  by  a  pair  of  plane  waves 
counterpropagating  in  the  ;  direction.  The  direction  of 
linear  polarization  of  each  field  is  shown  in  Fig.  1.  The 
total  electric-field  amplitude  is  given  by 

EU\r)  =  m(Ee-Ut+ik:  +c.c.)-fc2(Ff"“|-'K'  +c.c.). 


e-y2  e+y2 


e_3/2  e_y2  e+i/2  n+y2 


FIG.  2.  Atomic  level  scheme  and  Clebsch-t  Jordan  coeffi¬ 
cient*  for  fa)  ./  —  -t  '  K  ~  ~  and  (1)1  /.,  =  4  —  J.  —  4 
transitions. 


(2.1) 


where  the  real  amplitude  £'  does  not  vary  with  time. 
In  general,  the  unit  polarization  vectors  tq  and  v->  have 
different  directions,  so  that 


0  .  0 

tq  =  cos  —  er  +  sin  -oy , 
0  .  0 


e2  =  cos-cr  -  Sill  -ey 


(2.2) 

(2.3) 


where  the  angle  0  between  the  polarization  vectors  is  re¬ 
stricted  to  0  <  0  <  rr/2. 

The  laser  fields  interact  with  an  ensemble  of  atoms 
whose  ground  states  ij  are  characterized  by  total  angular 
momentum  J 9  =  1/2.  The  ground  states  are  linked  by 
(In'  nearly  resonant  laser  fields  to  excited  states  <  having 
total  angular  momentum  J,  —  1/2  or  3/2  [see  Figs.  2(a) 
and  2(b),  respectively].  Both  cases  are  similar  in  many 
respects.  To  be  specific,  in  the  qualitative  discussion  we 


FIG.  I.  One-dimensional  field  configuration  I  wo  line.ttlv 
polarized  fields  con nterpropagale  in  the  ;  direilion  with  an 
angle  H  between  their  polarization  vectors 


consider  the  case  Je  =  1/2. 

It  is  convenient  to  represent  the  field  (2.1)  as  a  sum  of 
two  circularly  polarized  fields  with  polarizations  such 
that 


<-±  = 


T«r  “ 


V~2 


and  amplitudes  L'±  given  by 
E±  —  T ''/2-t  cos  (k-;  ±  —  )  i 


(2.1) 


This  system  can  be  regarded  as  consisting  of  two  two- 
level  subsystems:  ;;  —  1/2  —  1 1/2  and  ;/ 1/2  —  <  —  1/2. 
driven  by  the  standing  waves  given  by  Fq.  (2.5)  and 
linked  to  each  other  by  the  transitions  involving  emission 
of  spontaneous  photons. 

The  force  exerted  hv  the  fields  on  the  atom  is  given  hv 


where  7/af  is  the  atom-field  coupling  given  by  Eq.  ( A 5 ) 
of  the  Appendix,  and  where  (  )  indicates  a  quantum  me¬ 
chanical  average.  In  the  weak-fieid  limit  defined  by  (1.1). 
one  can  adiabaticallv  eliminate  atomic  density-matrix  el¬ 
ements  involving  the  excited  state  and  obtain  the  force 
in  the  form 


f  —  c <  x !  -fi  /*,nr, 
where 


(2.7) 


/.V. 


ill  u 

i/r 


2  h<Uf- 
3(('1-  +  r-/ 1) 


cos  fi  sm(2f; ) 


(2.M 


and 
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/'’ml  =  /’+/'+  +  F-P- 
(P  +  P/4) 


sin  0  cos(2L-)Il . 


Tin*  quantities  appearing  in  Eqs.  (2.8)  anil  (2.9)  are  de¬ 
fined  as  follows:  the  Rabi  frequency  /  has  the  form 
/  =  /h ;  Vo  is  an  effective  potential  given  by 


/•+  and  /•'_  are  equal  to 

,,  =_F  '2hkfrf- 

+  3(6-  +  r’/'i) 


sin  0  cos(2ftc);  (2.11) 


p+  (p_  )  is  the  population  of  the  w j  —  1/2  ( w j  =  — 1/2) 
ground-state  magnetic  sublevel:  and  fi  is  the  population 
difference 

11  =p+-f>-.  (2.12) 

Consistent  with  the  weak-lield  limit  used  to  derive 
Eq.(2.7),  the  ground-state  sublevel  population  is  con¬ 
served,  i.e., 

(>+ +  P-  =  1-  (2.13) 

Equat  ion  (2.7)  is  valid  only  in  t  he  limit  of  relatively  small 
velocities, 

|/t|  <  max( E.  |A) ) ,  (2. I  f) 

which  is  the  only  limit  considered  in  this  paper. 

According  to  Eq.  (2.7),  the  force  consists  of  two  com¬ 
ponents,  an  “external"  part  /\.x,  which  is  independent 
of  the  internal  state  of  the  atom,  and  an  "internal"  part 
/•’mi,  which  depends  on  the  population  difference  II  oft  In- 
ground-state  sublevels.  The  internal  part  can  be  written 
in  terms  of  a  force  /-±  acting  on  each  sublevel  popula¬ 
tion.  On  averaging  the  force  over  an  optical  wavelength, 
the  contribution  of  the  external  part  vanishes,  while  that 
of  the  internal  part,  depends  on  the  Fourier  component 
of  If,  which  varies  as  cos(2Lc).  Thus  a  determination  of 
the  spatially  averaged  force  reduces  to  a  determination 

°N,(-~). 

The  population  difference  II  is  driven  by  t he  differential 
rati-  (!{+  —  li_)  at  which  atoms  are  optically  pumped  to 
the  +  and  —  ground-state  sublevels,  where,  as  shown 
explicitly  in  Sec.  III.  Ii±  are  given  by 


«±  =  7 


cos-(/-c  ±  ]■,()), 


(215) 


9(6,J  +  r-’/i) 


9(8-  +  l'7-l) 


and  f±  =  K±\tjh.  The  population  difference  relaxes  back 
towards  equilibrium  at  the  rate  ( I{+  -f  //_  ).  Fmler  the 
wenk-fiold  condition  (l.l).  in  steady  state  the  population 
difference  evolves  as 


i/II 

''Ih  ~(  i" 


V  )-(//*  I-  /«*_  )II . 


The  driving  term 

//+  —  //_  —  —  t~1  sin  </sin (m2kz)  (2.18) 

coincides  with  the  difference  of  light  shifts  of  ground-state 
sublevels.  It  is  a  periodic  function  of  atomic  position  and 
is  small  in  the  vicinity  of  nodes  (kz  «  rr / 2 -P  rr ? i )  and  anti- 
nodes  ( kz  ss  mi)  of  the  field  (the  positions  of  "nodes 
and  “antinodes”  to  which  we  refer  are  those  which  would 
appear  for  parallel-polarized  fields).  The  driving  term 
vanishes  for  0  =  0,  implying  that  suh-I)oppler-limit  cool¬ 
ing  does  not  occur  for  parallel  polarizations  of  the  fields. 
For  If  /  0  the  driving  term  increases  with  increasing  // 
and  reaches  its  maximum  value  for  laser  beams  having 
orthogonal  polarization  (0  =  -x). 

The  local  optical  relaxation  rate  /?|oc  for  the  popula¬ 
tion  difference  of  the  sublevels  is  equal  to 

/floe  =  «+  +  R.  =  T~ 1  [  1  +  cos  0  cos( 2/  ;  )].  (2.19) 

In  the  case  of  orthogonal  polarization  0  =  $ .  this  rate 
does  not  depend  on  position  c  However,  for  small  angles 
0  <C  1,  /fioc  is  strongly  position  dependent,  and  is  very 
small  near  the  nodes  of  the  field  (that  is.  relaxation  of 
the  ground-state  sublevel  population  difference  is  slow  in 
these  regions).  Although  the  driving  term  (2.18)  is  also 
small  in  t  he  vicinity  of  t  In-  nodes,  t  he  competition  of  t  hose 
two  small  quantities  leads  to  a  considerable  gradient  of 
the  population  difference.  It  is  the  rapid  variation  of  /f|,„- 
near  the  field  nodes  that  may  lead  to  qualitatively  new 
features  in  laser  cooling. 

To  understand  tin-  qualitative  features  of  tin-  spatially 
averaged  force,  we  firs!  consider  atoms  having  i  -  0.  Al¬ 
though  these  atoms  do  not  experience  the  averaged  force, 
it  follows  from  Eqs.  (2.7)  and  (2.17)  that,  for  sufficiently 
slow  atoms,  the  spatially  averaged  force  can  be  obtained 
in  terms  of  the  gradient  of  the  population  difference  of 
r  =  0  atoms.  From  Eq.  (2.17)  one  finds  this  population 
difference  to  be  given  by 


n(j-,c  =  o)  =  - 


sin  0  sin  2 kz 
1  +  cos  0  cos  2 kz 


For  0  =  x/2.  Il(r.r  =  0)  =  —  siu(2L:).  Defining  a  char¬ 
acteristic  scale  ;M.  over  which  II  varies  as 

-s<-  55  1 1  /( c/H  / </-  )-  (  2  2  I ) 

one  sees  that  cs,.  ~  k~]  when  0  =  x/2  (see  curve  I  in  Fig. 
3).  However,  for  small  angles  f)  «  I .  a  new  length  scab 
is  introduced  into  the  problem.  Far  from  I  lie  field  nodes 
( \kz  -  (j (  1  +  2n ) |  ^  I) ,  n  =  0.  ±  I .  .  )  l  In-  populat  ion 
difference  is  small: 

life,  r  =  0)  =  —0  tan  kz,  (2.22) 

while,  in  the  vicinity  of  the  nodes,  ll(r.r  =  0)  takes  the 
form 

where  <  =/':  —  (  I  +  2n )  x/2  One  can  see  from  Fq  (2  23) 
and  curve  2  in  Fig.  3  that,  in  contrast  lothed  =  x/2  case, 
for  small  0,  near  the  nodes  the  population  difference  is  of 
order  unity  and  varies  very  rapidly  on  a  scale  ;S1  ~  2L.  i 
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FIG.  3.  The  population  difference  1J(-,  u  =  0)  of  the 
ground-state  sublevels  as  a  function  of  A-  for  t)  =  -f  and 
0  —  0.2.  The  value  k:  --  -f  corresponds  to  a  field  node. 


Next  we  consider  an  atom  moving  in  the  ;  direction 
with  a  very  small  velocity  v,  such  that  11  does  not  vary 
considerably  during  an  optical  pumping  period,  that  is, 

«'  <  /?ioc~%c-  (2.24) 

In  this  limit  the  velocity-dependent,  part  11,  of  the  pop¬ 
ulation  difference  grating  which  gives  rise  to  a  non¬ 
vanishing  averaged  force  can  be  obtained  from  Eq.  (2.17) 
by  a  standard  one-step  iteration,  and  is  given  by 

_  ,-i  dll(?,  r  =  0)  _  2Aer!,(eos0  -f-  cos'2kz) 

v  ~  V  loc  dz  ( 1  +  cos  0  cos  2Ar)3 

(2.25) 

Atoms  having  velocities  r  satisfying  ine(|iiality  (2.24)  give 
rise  to  a  spatially  averaged  force  that  varies  linearly  with 
v  (so-called  “capture  range”  of  the  force  [13]).  One  can 
see  from  Eq.  (2.25)  that  a  decrease  in  the  optical  relax¬ 
ation  rate  and  an  increase  of  the  population  difference 
gradient  lead  to  larger  II,  (see  Fig.  1)  and  might  thus 


FIG.  -I.  Velocity-dependent  part,  of  the  population  differ¬ 
ence  1 1  (  -.  r)  that  contributes  to  the  spatially  averaged  force, 
as  a  function  of  k:  for  several  values  of  H  and  o  =  krrv.  (I) 
H  =  -s.o  =  l.ll  X  Itr';  ( J)  H  =  (1.2.  O  =  1.0  X  10-’;  (  t) 
0  =  j,n  =  r,  X  It)—';  (  I)  0  =  0.2.  n  =  a  x  It)-'. 


lead  to  stronger  averaged  force.  Owing  to  the  depen¬ 
dence  of  R\oc  on  0 ,  both  the  capture  range  defined  by 
(2.24)  and  the  magnitude  of  the  force  depend  strongly 
on  0. 

For  the  fields  with  orthogonal  polarization  (0  =  it/2). 
one  has  R^c  —  rp,  and  E'q.  (2.24)  takes  the  form 

o  =  kvTp  <  1.  (2.20) 


In  this  limit  one  obtains  the  well-known  result  [1]  for  the 
friction  force 


25  A-/1' 6a 
3(62  +  r-’/4)' 


(2.27) 


For  small  0,  both  the  local  optical  relaxation  rate 
and  the  spatial  scale  over  which  the  population  differ¬ 
ence  varies  decrease  near  the  field  nodes,  so  that  /?|oc  ~ 
Tpl0~/'2  and  ;sc  ~  0/2.  It  follows  from  Hqs.(2.19)  and 
(2.24)  that  the  capture  range  is  given  by 

a  C  03  (2.2b) 


and  is  much  smaller  than  that  for  0  =  r; .  In  this  capture 
range,  however,  the  st  rong  gradient  of  II  and  long  optical- 
pumping  time  in  the  vicinity  of  the  field  nodes  give  rise 
to  a  spatially  averaged  force 


25  A  /'6a 
30(6-’+  I'-’/l)  ' 


(2.29) 


which  is  larger  than  that  for  0  =  For  small  0  r.nd  t 
the  entire  contribution  to  the  force  originates  when  an 
atom  moves  near  the  field  nodes. 

The  iterative  approach  used  to  arrive  at  Eq  (2.25)  is 
only  applicable  when  condition  (2.24)  is  satisfied.  In  the 
case  of  orthogonal  field  polarization,  condition  (2.24)  is 
satisfied  provided  that  a  -  kirr  <§;  I;  however,  for  0  -C  1 
this  condition  is  satisfied  only  for  a  =  Apr,,  <C  03  Thus 
for  0  <  1  there  exists  a  range  of  velocities 

l»n»  03,  (2.30) 


for  which  the  iterative  approach  fails.  These  velocities 
are  high  enough  to  prevent  the  population  difference  of 
the  ground-state  sublevels  to  undergo  considerable  re¬ 
laxation  during  the  passage  through  the  field  node,  but 
low  enough  for  this  process  to  occur  at  a  spatial  scale 
smaller  than  the  optical  wavelength.  For  this  velocity 
range,  a  theoretical  approach  based  on  a  Fourier  series 
expansion  for  II  is  not  of  much  practical  use  either,  since 
the  number  of  terms  involved  is  prohibitively  large.  New 
analytical  techniques  are  needed,  which  are  described  in 
tin'  subsequent  sections. 

Nevertheless,  it  is  possible  to  obtain  a  qualitative  pic¬ 
ture  of  the  velocity  dependence  for  the  force  in  the  range 
(2.39)  by  introducing  the  concept  of  an  effective  optical 
relaxat  ion  rate  /f,  ir  denned  by 

/fto,(0*fl'=  1.  (2  31) 

I 
l 

where  I  —  (z  —  :o)/p.  and  R |,„  is  given  by  lap  (2  19) 
Equation  (2.31)  simply  states  that  the  probability  of  op¬ 
tical  pumping  in  the  time  interval  (/  -  .  I )  is  equal 
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to  unity.  For  the  range  of  atomic  velocities  satisfying  in¬ 
equalities  (2.30),  in  the  vicinity  of  the  nodes,  integration 
of  Eq.  (2.31)  leads  to  the  effective  optical  relaxation  rate 

fleff  ~  r - 1  [c-  +  (jjo )2/3]  «  T~ 1 ,  (2.32) 


which  does  not  depend  on  0  and  increases  with  increasing 
velocity.  The  spatial  scale  of  the  population  difference  is 
now  given  by  csc  =  v R~j  ~  al^3/k  (see  curve  4  in  Fig. 
4).  The  largest  contribution  to  the  spatially  averaged 
force  is  produced  in  a  single  optical-pumping  cycle  when 
an  atom  passes  the  field  node.  In  the  range  of  velocities 
(2.30)  the  averaged  force  is  given  by 


tik  f~6  sin2  0 
3v/3(^2  +  r2/4) 


(2.33) 


and  is  independent  of  v. 

A  graph  of  F  vs  0  for  small  a  -C  1  is  shown  in  Fig.  5. 
The  force  first  increases  with  increasing  0,  and  reaches 
its  maximum  value  for  a  small  angle  0  ~  o1/3.  Then 
it  decreases,  approaching  the  value  F(0  —  t/2,  v)  as  0 
tends  to  7t/2. 

We  now  turn  our  attention  to  the  momentum  diffusion 
coefficient  D ,  which  consists  of  two  parts: 


D=D,p  +  Dind,  (2.34) 

The  first  contribution  Dsp  is  related  to  the  fluctuations 
concerned  with  emission  of  spontaneous  photons.  This 
contribution  is  quite  similar  to  the  momentum  diffusion 
coefficient  in  a  two-level  system.  For  |Arr|  T,  it  de¬ 
pends  only  slightly  on  angle  and  velocity,  since  the  spe¬ 
cific  features  of  sub-Doppler-limit  cooling  and  mult  ilevel 
dynamics  do  not  play  a  significant  role  in  the  diffusion 
processes  brought  about  by  these  fluctuations.  The  av¬ 
eraged  contribution  Dsv  can  be  roughly  estimated  by 


(2.35) 


The  second  contribution  D,,„i  is  produced  by  the  lluc- 
tuations  of  the  stimulated  force  exerted  on  the  atom,  and 
is  given  by 


FIG.  5.  The  spatially  averaged  laser  cooling  force  as  a 
function  of  angle  6  for  small  n  =  krrp  =  1 .11  x  I'l-’. 


[(F(t)F(t  -  t))  -  (F(t))(F(t  -  r))}dr . 


This  contribution  is  directly  related  to  the  multilevel  dy¬ 
namics,  since  it  is  brought  about  by  fluctuations  in  the 
time  periods  an  atom  stays  in  sublevels  mj  —  ±1/2  of 
the  ground  state.  For  relatively  small  velocities  satisfying 
the  condition  q  Cl,  the  correlation  time  of  the  forces 
in  Eq.  (2.36)  is  given  by  R~ft f,  and  D,nd  is  approximately 
equal  to 


where  F+  is  given  by  Eq.  (2.11). 

For  9  <C  1,  the  force  E+  is  9  times  smaller  than  that 
for  0  =  f.  On  the  other  hand,  in  the  vicinity  of  the  field 
nodes,  <C  r  ~ 1 .  The  overall  effect  is  a  decrease  of  the 
averaged  diffusion  coefficient  with  decreasing  6.  given  by 

96-  - 

Ojnd((?,  r)  - — pY DSp  (2.38) 

for  very  slow  atoms  satisfying  Eq.  (2.28),  and 


And(fU') 


o>/3r2' 


for  intermediate  velocities  (2.30).  Thus,  for  small  ft.  the 
induced  momentum  diffusion  coefficient  D,„ rj(0  <C  I  t  ) 
is  much  smaller  than  Din d{0  =  tt/2,  t  ).  It  also  exhibits 
significant  dependence  on  velocity  v  in  the  range  O'  <c 
(>Cl.  The  consequences  of  the  9  dependences  of  F  and 
D  are  discussed  below. 


III.  EQUATIONS  FOR  THE  GROUND-STATE 
SUBLEVEL  POPULATIONS 


In  this  section,  using  the  Wigner  representation  [14], 
we  derive  the  equations  for  the  populations  of  the  ground- 
state  sublevels.  We  consider  atoms  having  ground-  and 
excited-state  angular  momenta  J y  ^  and  Jr  =  i  or  f 
Using  the  general  equations  (see,  for  instance,  [14,4])  for 
the  density  matrix  p(:,p).  where  p  =  Mv  is  the  atomic 
moment  um,  and  c  is  the  center-of-mass  posit  ion,  and  tak¬ 
ing  into  account,  the  weak-field  condition  (1.1),  one  can 
adiabatically  eliminate  excited-state  matrix  elements  and 
obtain  a  closed  set  of  equations  for  the  density-matrix  el¬ 
ements  of  the  ground  state  [1,4].  The  derivation  of  (In¬ 
equations  is  given  in  the  Appendix  and  is  quite  similar 
for  the  two  values  of  Assuming  that  the  atomic  mo¬ 
mentum  is  much  larger  than  that  of  the  photon,  that  is 

p  >  hk,  (31) 

one  finds  the  following  equations  for  the  ground-state 
sublevel  populations: 
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dtP±  9(«2  +  P/4) 


for  Je  =  k  and 


o  (9  ,  \  /  ,  0  \  3 hkb  . 

-  cos-  /  T-  +  k;  J  p±  +  cos'  I  ±-  -f-  k ;  J  - —  si n(^fO  +  2b:) 

+  ~~  f[M  -  11  cos(T0  +  2 kzj}^-f±-  +  [(i  +  cos(±0  +  2b:)]^£ 
20  V  (//>-  op~  / 


Jtp±  =  WTWJa)  ["cos"  It2  +  k:)  p±  +  cosJ  +  k:)  ^  -  -^f [sin(^  +  2A';)  +  :is,n(±* +  ‘2K')]17 

/  ~  lc‘"~  (  *V2  ^2  \ 

+  ■' —  (  [70  -  13cos(:y0  +  21c;)  -  27cos(±0  +  2 k;)]-:-P^  +  2[6  +  cos( ±0  +  2k;))—~~-  ) 
40  V  op-  tip-  / 


3 hkb.  ,dp+ 

p? - — —  [sinful?  +  2  icz)  +  3sin(±0  +  2k;)}— — 


for  Je  =  f,  where  d  /dt  =  d  /dt  +  vd  /d:.  Equations 
(3.2)  and  (3.3)  provide  a  complete  description  of  the  time 
evolution  of  the  atomic  distribution  of  ground-state  sub¬ 
levels  as  a  function  of  atomic  center-of-mass  position  and 
momentum,  in  the  limit  of  weak  field  (1.1)  and  small  pho¬ 
ton  momentum  (3.1).  The  first  two  terms  on  the  right- 
hand  side  of  the  equations  describe  the  spatially  tie  pen¬ 
dent  population  exchange  between  the  sublevels  owing  to 
the  optical  pumping.  The  term  with  d/dp  can  be  inter¬ 
preted  as  a  gradient  of  an  effective  potential  determined 
by  the  spatially  dependent  ac  Stark  shifts  of  the  levels. 


for  J,  =  4-  and 


-[cos(T0  +  2k;)  +  3cos(±0  +  2b:)] 


1 2  ( A  -  -f  r-‘/4) 


for  ,7 e  =  fj .  Finally,  t  he  terms  containing  the  second 
derivatives  with  respect  to  p  correspond  to  the  momen¬ 
tum  diffusion  process  caused  by  emission  of  spontaneous 
photons. 

To  find  the  atomic  momentum  and  position  distribu¬ 
tions  it  is  useful  to  introduce  the  probability  density 

t’U.P)  =  f>+  +  P-  (T(>) 

for  an  atom  to  have  momentum  ;>and  position  :.  and  the 
population  difference  density 

II (;■}>)  =  P+ ~  P--  (d.7) 

Then,  using  Fq.  (3.2)  with  .7,  =  4,.  one  arrives  at 


cos 0 sin  21c z— —  sin  0  c os  2k; - 

Up  (ip 


(2  —  Co;  rnv  2k;)- 


-  sill  0  sin  2b;- 


—  n  =  — -  —V’  sin  0  sin  2k;  —  (1  +  cos  0  cos  2b;)fl 

dt  tp 

Zhkb  (  dil>  .  dn 

H - sm  0  cos  2b:— - cos  Osin  2b;— — 

r  V  dp  dp 


frk2  (  „  „  d-ri 

H - 1(2  —  3  cos  0  cos  2k;)  -— 

5  \  dp 2 


—  3  sin  0  sin  2b;- 


wliere  the  “optical  pumping"  time  rp  is  delined  by  Eq. 
(2.1(i).  Similar  equations  can  be  obtained  for  .7,  =  4 
using  Eq.  (3.3). 

In  general,  Eqs.  (3.2)  and  (3.3)  or  (3.8)  and  (3.9)  can¬ 
not  be  further  simplified.  However,  if  the  atomic  kinetic 
energy  Ek  —  p~/2M  is  much  larger  than  the  effective 
potentials  S±,  that  is  if 

jg»S»£rW  (:U0) 

the  atoms  move  almost  freely  and  localization  effects  are 
weak.  This  implies  that  the  time  scale  rp  that  character¬ 
izes  the  time  variation  of  r(;./>)  satisfies  the  condition 

Tpkv  1  (3.11) 

As  a  consequence,  the  probability  v[;,p)  —  p+  + /'_  de¬ 
pends  only  slightly  on  :.  Below  we  neglected  this  de¬ 
pendence  [15]  The  problem  is  then  reduced  to  a  deter¬ 
mination  of  the  atomic  momentum  distribution  .(/>)  = 
V'f-./')-  To  obtain  i/'(p)  as  a  solution  of  Eq.  (3.8),  one 
must  know  the  population  difference  density  Il(:./>)  The 
steady-state  distribution  for  the  population  difference 
II ss  is  established  on  a  time  scale  H~J  given  In  Eq 
(2.31).  Although,  as  has  been  discussed  in  Sec  II.  /i,u 
may  vary  significant ly  with  atomic  position  am  I  velocity. 
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it  can  be  shown  ( see  below)  that  /p  <C  r/.;.  In  detormiii-  terms  containing  OU/dp.  and  d~y</dp~ and  cl- II /Op-  in 
ing  if'(p)  from  Eq.  (3.S),  one  is  then  justified  in  replacing  Eq.  (3.9)  [owing  to  conditions  (3.1)  and  (3.10),  respec- 
II  m  that,  equation  by  its  steady-state  value  llss.  Tak-  lively],  one  finds  the  steady-state  solution  of  Eq.  (3.9)  in 
ing  into  account  the  condition  E\.  ;§>  S±  and  neglecting  the  form 


llss(-.p)  —  f  t lz'  (- sindsin(2/;r  -  z')iiiz,p)  +  — -7-  sin  0  eos(2kz  -  j— 1 -  * 

2a  Ja  V  1  (>V 


x  exp  — (2a )“ 


+  cos  0  eos(2kz  —  ;")](/;' 


a  =  lcvTp.  (3.13) 

Substituting  H ss ( -  •  p)  in  the  form  (3.12)  into  Eq.  (3.8)  and  neglecting  derivatives  of  order  higher  than  2.  one  arrives 
at  a  Fokker-Planck  equation  of  the  form 

T,v  =  4-1- F  +  |-(0,PV)]  ■  (3.14) 

ilt  op  dp  dp 

where  the  force  F  is  given  by 

F  =  3(p‘+r'-/4)  [cos£,sill2i”  +  -sin-  cos 2Ar 

x  J  clz'  sin(2&;  —  ;')  exp  (2a)-1  J  [1  +  cost?  cos(2kz  —  r")]d;"^  ,  (3.15) 

tiie  sub-DoppIer-limit  “stimulated”  contribution  Dm(\  to  the  momentum  diffusion  produced  by  the  fluctuations  of  the 
instantaneous  dipole  forces  by 

Dlnd  =  ~ — ^  cos2kz  J  dz'  cos(2kz  -  z')exp  ^-(2o)_1  J  [1  +  cos0  cos{2kz  -  c")]r/r"^  ,  (3.16) 
and  the  “spontaneous”  contribution  Dsp  to  the  momentum  diffusion  associated  with  emission  of  spontaneous  photons 

by 


,,  h  n  n  sin  Osinilrr 

l)s p  = -  2  —  cos  0  cos  2k  z  + - - - 

2  rp  2a 


—  J  dz' s\n(2kz  —  c')exp  ^-(2a)  1  J  [1  -I-  cos0cos(2kz  —  c")]c/c"^ 


Under  assumption  (3.10)  considerable  variation  in  atomic  kinetic  energy  occurs  on  a  time  scale  larger  than  (Ar)-1 . 
In  other  words,  condition  (3. 1 1)  is  satisfied.  In  this  limit,  the  force  and  diffusion  averaged  over  a  wavelength  determine 
the  time  evolution  of  the  distribution  ip.  Averaging  Eqs.  (3.15)-(3.17)  over  a  wavelength,  and  substituting  the  results 
into  Eq.  (3.14),  one  can  arrive  at 

»  =  +  +  !;(«,„«].  <*»> 
where  the  spatially  averaged  force  /  is  given  by 

r  _ _ — r  1  (**0  _  \ 


/••  =  -- 


ill -0  p  .  _r,  feast)  \  , 

-)Jo  l^-—su>(nr))dr. 


3( A-  +  P-/4)  ,/(l  '  v 

and  the  averaged  diffusion  coefficients  an 

„  v  r,  ,1  :  r 


D„,,  = 


2h  k-fJV  (  1  r  _  r  frost)  \ 

=  W  +r-7D  V  +  V"r0L  [«-(2-)/..(— s„,(-))- 

U~k2ti'f- sill"  0  f'  \  feast)  \  feast) 

"WTrwi  ''  . a) +  /=(—■ 


(,ir0  —  /,  ( - 'in(itr)]  ,/rl  ,  (II ’.'ll) 


111(0  r)  Jr 
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where  I0  and  72  are  modified  Bessel  functions  [16]. 

For  Je  =  one  can  also  obtain  equations  similar  to  Eqs.  (3.14)  and  (3.18).  For  example,  the  force  F  is  given  by 


F 


hk6f 2 


3(  <82  +  T2/4) 


1 2  cos  0  sin  2 kz  —  (4a)  1  sin2  0  cos  '2k: 


x  f  dz'  sin{'2kz  —  c')exp  (  —(4a) 
Jo 


-  / 1. 


+  cosO  cos(2kz  —  z")]dz" 


and  the  averaged  force  F  by 
hk  f  ~6  sin2  0 


F  = 


sin(2ar)  J  dr  . 


6(^2  +  r2/4)  J0  '■  V  2« 

For  the  sub-Doppler-limit  stimulated  part  of  the  momentum  diffusion  coefficients  one  finds 


D,"d  =  2^  cos(20-.-')exp(-(4a)->|  [1  + 


>s0cos(2/fcc  —  z")\dz"  dz',  (3.24) 


(3.22) 


(3.23) 


hrk-Pp  sm-e  f00  _T 


And  2r^2  +  p/4)  l  C 


,  .  cos  0  .  \  ( cos  0  . 

cos(4ar)/0  (  — — sin(2ar)  1  +  /j  I  — — sm(2ar) 


dr. 


while  the  spontaneous  contribution  Dsp  to  the  momentum  diffusion  is  given  by 


A 


sp 


trk2 

40r,, 


41  -  19  cos  0  cos  21:;:  — 


2  sin2  6  sin  2  kz 


J  dz' sin(2fcz  —  .s')  exp  /  —(4a)  1  J  [1  +  cos(?cos(2£z  —  z")]dz' 


and 


f  c  T  cos(4or)/u  (  C°^  sinf. 

do  L  V  'd<> 


2ar)  ]-/■’(  sin(2or) 
2a 


(3.25) 


(3.26) 


dr\  .  (3.27) 


IV.  DISCUSSION 

In  this  section,  if  not  stated  otherwise,  we  consider 

Je  =  i 


A.  Force  and  diffusion 

For  0  =  ~  one  can  use  Eqs.  (3.19)  (3.21)  to  obtain 
expressions  for  the  force  and  momentum  diffusion  coeffi¬ 
cients  in  a  closed  form  for  arbitrary  values  of  the  dimen¬ 
sionless  velocity  a  =  krrp.  For  the  averaged  forte  one 
recovers  the  well-known  result  [1] 


'(*-!)- 


2 hk f-fin 


3(A+  r'2/4)(l  +4o3)' 


The  diffusion  coefficients  are  given  by 

frk-f>-f- 


AikI  (0  =  J) 


n(f>~  +  rL’/I)(  I  +  la2) 


4.1) 


;i.2) 


and 


D. 


sp 


(•=§) 


2trk-f-r 

m2  +  r2/4) 


1  + 


1 


1(1  +  la'-’) 


(4.3) 


and  their  qualitative  behavior  coincides  with  that  of  dif¬ 
fusion  coefficients  obtained  in  Ref.  [4] 

For  an  arbitrary  angle  0.  it  is  possible  to  obtain  an¬ 
alytical  expressions  for  the  averaged  force  and  diffusion 
coefficients  in  the  limiting  cast's  of  small  and  large  veloc¬ 
ities.  Numerical  solutions  for  the  force  anti  the  diffusion 
coefficients  as  a  funct  ion  of  a  are  presented  in  Figs  ti  and 
7,  respectively. 

1.  Low  velocities  at  <C  1 


In  this  limit,  for  verv  low  velocities 


03 

o«- 


(4.4) 


one  can  take  the  integrals  in  Eqs  (3.19)  (3.21)  to  obtain 
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f 


F  —  — 


'2hk  f26a 

3  sin  0(62  +  r'-/4) 


2 hkf-6a 
'3(?(62+  r2/4) 


3/i  k-bv 


D  —  Dsl  +  Ds0  — 


h2k2f~  r 

9(^2  +  r2/4) 

/  sin2f? 


x  12+tt^  + 


1862  sin  0  \ 

r2  i  +  sin  o ) 


The  force  increases  and  the  diffusion  coefficient  decreases 
gradually  with  decreasing  9. 


D  2h2k2}2Y  /  9  S29\ 

9(^2  +  r2/4)  v  +  r2  ; 


The  force  increases  as  0~l  and  the  diffusion  coefficient 
decreases  linearly  with  decreasing  0. 

To  determine  the  source  of  such  a  strong  dependence 
of  the  force  and  diffusion  on  0,  one  can  examine  the  ex- 


In  particular,  for  small  angles  6  <C  1  and  very  small  pressions  (3. 15)— (3. 17)  for  the  nonaveraged  force  and  dif- 


velocities  satisfying  (4.4),  one  finds 


fusion.  If  condition  (4.4)  is  satisfied,  one  has 


„  2 hk6f-  A  .  ,  .  sin  4k; 

F  sa - - - — —  cos  0  sm  2k;  +  sin"  0\  - - - - — - 

3(62  +  1  2/4)  \2(  1  +  cost?  cos  2k:) 


2(.\  cos  2 kz( cos  2 k:  4-  cos  0 ) 
( 1  +  cos  0  cos  2 k; ):! 


(b) 


FKi.  0.  Spatially  averaged  laser  cooling  force  as  a  function 
of  o  =  kt’Tp.  The  different  curves  correspond  to  dillerenl  an¬ 
gles  0  between  the  fields'  polarization  directions.  The  atomic 
transition  occurs  between  two  levels,  each  having  total  angu¬ 
lar  momentum  ./  =  and  0  =  51'.  The  force  for  very  low 
velocities  is  shown  in  more  detail  in  (b). 


FIG.  7.  The  spatially  averaged  momentum  diffusion  co>  f- 
ficicnt  as  a  function  of  o  =  kvrp.  The  different  curves  cor¬ 
respond  to  different  angles  0  between  the  fields  polarization 
directions.  The  atomic  transition  occurs  between  two  levels, 
each  having  total  angular  momentum  ./  =  4.  and  0  =  51' 
The  momentum  diffusion  coefficient  for  very  low  velocities  is 
shown  in  more  detail  in  (b). 
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For  0  -C  1,  far  from  the  nodes  of  the  standing  wave  [|&-  — 
^(1  +  2?i)l  ~  1]  the  velocity-dependent  part  of  the  force 
(4.9)  is  of  order  02hk26v/T  and  is  much  smaller  than 
that  of  the  force  obtained  for  0  =  t/2.  However,  in  the 
vicinity  of  the  nodes,  where 

C=:jfes-|(l  +  2n)«l,  (-1.10) 

the  force  takes  the  form 

„  2  likSf2  /2f(0-/4  —  f2)  02cv(02/4 -<2)\ 

'  %  +  r-'/d)  V  (<-  +  0-74)  27-  +  0-7-ip  J 

(in) 

The  first  term  in  Fq.  (4.11)  is  velocity  independent  and 
corresponds  to  the  gradient  of  the  effective  potential  G,,ir 
given  by 

i'-  =  3F^(<5-Ti’''',  +  i’T))-  («■»> 

The  second,  velocity-dependent,  term  can  In-  seen  in  Fig. 
-I  (curve  2).  In  the  very  vicinity  of  the  node  <  0.  it 

represents  a  very  strong  frictional  force 


Ev 


'Mhk-6v 

fl-T 


(114) 


while  for  |c(  >  0/2  it  changes  sign  and  leads  to  heating. 
However,  after  averaging  over  the  position  around  the 
node  the  combined  effect  is  still  that  of  cooling  given  by 
Eq.  (4.7).  Thus  one  can  see  that  for  small  0  a  strong  fric¬ 
tional  force  is  pioducetl  in  vicinities  close  to  the  standing 
wave  nodes. 

'The  non  averaged  sub-DoppIer-Iimit  stimulated  contri¬ 
bution  /};„<)  to  the  momentum  diffusion  is  proportional 
to  O'-  far  from  the  nodes  and  is  much  smaller  than  that 
in  the  case  0  =  t/2.  However,  near  the  field  nodes  [i.e., 
when  condition  (4.10)  is  satisfied],  Dm,\  is  given  by 


1  /Uik(-0\-  1 

D"ut  ~  2r,,  V  F  J  c-  +  02 74 


(4.14) 


For  (  <  0 ,  l)„i,\(0  1)  is  of  the  order  of  D,w\{0  = 

t/2)  and  provides  the  major  contribution  to  the  averaged 
diffusion  coefficient.  (4.16). 

For  the  intermediate  range  of  velocit  ies 


<P 

1  >  o  >  — 

6 

the  averaged  force  takes  the  form 

/■  -  2hk 

~  4v/4(/>2  +  1-/4)' 


(4.15) 


(4.16) 


Thus,  for  <«<l  tin*  force  (-1.16)  does  not  depend 
on  tin'  velocity  and  decreases  with  decreasing  0.  The 
diffusion  coefficient,  varies  as 


„  Irk2  r-Y  /  9.54620-\ 

1  ~  <)(fi2  +  f-/4)  (2+  r2«>/»  )  '  '  ‘  ° 

and  decreases  with  decreasing  0  and  increasing  velocity. 


2.  Large  velocity,  or  1 


In  the  high-velocity  limit,  the  force  takes  the  form 
0hkf26  sin2  0hfAV6  s'nr 


F  =  - 


6  a(6-  +  F2/4)  27  v(S2  +  r2/4)2  ’ 


(4.18) 


which  is  valid  for  any  angle  0.  The  force  (4.18)  decreases 
with  decreasing  0. 

The  diffusion  coefficient  is  given  by 


D  = 


2trkrf-Y 
9(P  +  F2/4) 


(4.19) 


and  decreases  with  decreasing  0  and  increasing  velocity. 

Results  for  the  entire  range  of  velocit  ies  are  shown  in 
Figs.  6  and  7.  The  velocity  dependence  of  the  force  in 
the  case  of  small  angles  0  differs  dramatically  from  that 
obtained  in  the  case  0  —  t/2.  For  a  given  o  the  force 
F(0  =  t/2)  is  smaller  than  F(t/2  >  0  >  o'^2).  An 
important  new  feature  of  the  diffusion  coefficient  is  the 
strong  dependence  of  its  sub-Doppler-limit  part  Djnd  on 
angle  0.  For  0  =  t/2  and  |/>|  ;§>  F  this  part  is  the  domi¬ 
nant  one,  as  it  is  A~/I'2  times  larger  than  the  contribution 
Ds p  brought  about  by  the  emission  of  spontaneous  pho¬ 
tons.  However,  /And  decreases  with  decreasing  0 ,  while 
Ds p  does  not  vary  significantly  with  0.  As  a  result.,  the 
diffusion  coefficient  D  decreases  with  decreasing  0.  More¬ 
over,  for  a  range  of  angles  0  1.  one  finds 


Duu\  <c  Dsl> 


(4.20) 


even  for  small  velocities. 


B.  Equilibrium  momentum  distribution 


For  J,  =  ^  and  />  >  0,  the  equilibrium  momentum 
and  spatial  distribution  t!vq(::,/>)  characterized  by  mean 
kinetic  energy  E =  p2u/(2Af)  results  from  the  balance 
between  the  cooling  force  and  diffusive  heating.  Rigor¬ 
ously,  the  equilibrium  distribution  is  a  solution  of  Eqs 
(4.2)  and  (4.4);  however,  it  may  be  possible  to  approxi¬ 
mate  this  distribution  using  the  Fokker-Planck  equation 
(4.18)  with  a  diffusion  coefficient  D(0,p)  and  a  drift  term 
F{0,p).  The  resulting  distribution  t7q(/>)  is  not  necessar¬ 
ily  of  Gaussian  form.  Introducing  a  dimensionless  atomic 
momentum 


= 


_P_ 

t,k' 


(4.21) 


.ind  neglecting  the  dependence  of  Dsp  on  velocity,  one 
finds  the  distribution  tVq ( 7 )  to  have  the  form 


VVq(<7  =  VVq(0)e.\p 


(4.22) 


Calculating  the  mean  equilibrium  k  netic  energy  A1,-.,  as 


E =  En  :  /i'7y.q  (//)(///, 

J  -.X) 


(4.24) 


where  En  =  (hk)2 /(2M )  is  the  recoil  energy,  one  arrives 
at  the  results  presented  in  Figs.  8  and  9. 

’file  mean  kinetic  energy  as  a  function  of  dimensionless 
laser  field  intensity  /  defined  by 
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FIG.  8.  Mean  equilibrium  kinetic  energy  in  units  of 
h‘k2/2M  (M  is  atomic  mass)  as  a  function  of  a  dimensionless 
field  intensity  /  defined  by  J  =  |/|2A//(3h6Jt2).  The  different 
curves  correspond  to  different  angles  8.  The  atom-field  con¬ 
figuration  is  the  same  as  that  used  to  obtain  Figs,  ti  and  7 


7  =  W  =  0,.-  =  0)/2£„  =  Li|. 

is  shown  in  Fig.  8  for 


(‘121 ) 


x/-l>0>V/6.  (4  25) 

Assuming  that  |6|  3>  T,  and  using  expressions  for  /•'  and 
D  in  different  velocity  limits,  one  can  get  a  qualitative 
understanding  of  the  results  presented  in  this  figure.  For 
0  =  7r/2  one  arrives  at 


t/v  q(/?)  =  Vw(0)(l  +  2/?2/ r-)-"A  (4.26) 


and 


f‘Jeq  — 


Er 


r- 

I  -  6' 


(4.27) 


If  6  <  tr/2,  for  relatively  small  intensity  /  satisfying  the 
condition 


I  >  6  sin  2  6  , 


(4.28) 


one  has 


</-«,(/?)  =  0«,(O)[l  +  2/?7(/sin0)2] 


2l  — /  sin  J  C/4 


and 


/2sin2  0 


(4.29) 


(4.30) 


For  0  1  and  /  satisfying 


(&)'»>  Off- 

(4.31) 

one  arrives  at 

(4.32) 

and 

(4.33) 

For  even  higher  intensities 

(4.34) 

the  energy  distribution  takes  the  form 

V'Vq  =  <M0)exp(-(?2/2/ sin2  6), 

(4.35) 

and  the  mean  kinetic  energy  is  given  by 

Ee q  =  ErI  sin2  0. 

(4.36) 

One  can  see  from  Fig.  f  that  the  minimal  intensity  that 
can  still  ensure  an  equilibrium  energy  distribution  in¬ 
creases  with  0  as  6/sinJ0.  The  optimal  cooling  is  ob¬ 
tained  for  0  =  and  /  =  12.  However,  for  smaller  0  and 
higher  intensities,  the  mean  kinetic  energy  as  a  function 
of  /  increases  more  slowly  than  for  0  =  it/l,  and  in  some 
angle  range  can  be  a  decreasing  function  of  0.  The  de¬ 
pendence  of  Ee q  on  0  is  shown  in  Fig.  9.  For  relatively 
low  intensify  satis  ying 

(4.37) 

E„ q  is  given  by  Eqs.  (4.30)  and  (4.33)  and  is  almost  in¬ 
sensitive  to  the  angle  unless  it  is  very  small.  However, 
for  higher  intensity 


FIG.  9.  Mean  equilibrium  kinetic  energy  ill  units  of 
fi2k2 /2M  (M  is  atomic  mass)  as  a  function  of  8.  The  different 
curves  correspond  to  a  different  dimensionless  field  intensity  / 
defined  by  /  =  |/)2M/(3fi6A-2).  The  atom-field  configuration 
is  the  same  .is  that  used  to  obtain  Figs.  G  and  7.  and  angles  satisfying 


(4.38) 
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tt/2  >  0  » 


(1.39) 


the  equilibrium  energy  is  given  by  Eq.  (4.36)  and  varies 
as  sin'  0.  If  the  angle  gets  even  smaller,  that  is,  if 


»  0  »  0.5 


(4.40) 


/v0(|  is  given  by  Eq.  (4.33)  and  does  not  depend  on  0. 
Finally,  for  very  small  angles 


0.5  | 

f  r* 

\Sr>I 

one  has 

Ee q 

=  El 

1/8  /fl\l/2 

»  0  > 


(?) 


^7777 


8/2r2 


962(/02  -  4)(/02  -  6) 


(4.41) 


(1.42) 


energies  (4.44)  and  strong  localization.  For  Je  =  1/2  the 
potential  f/o  has  its  minima  at  the  field  nodes.  Thus  an 
atom  trapped  in  a  region 


near  the  field  node  may  experience  a  very  strong  fric¬ 
tional  force  and  may  be  cooled  to  the  energies  given  by 
Eq.(4.44).  However,  the  effective  potential  Uq  of  the 
Je  =  3/2  atom  has  its  minima  at  the  antinodes  of  the 
laser  standing  wave.  As  soon  as  an  atomic  kinetic  energy 
becomes  smaller  than  2(/o(z  =  0),  an  atom  is  trapped  in 
the  vicinity  of  the  antinode  where  both  the  force  and  the 
sub-Doppler-limit  diffusion  coefficient  are  proportional  to 
02  and  are  small.  As  a  result,  one  should  have  been  ex¬ 
pect  to  obtain  weaker  localization  (6z  ~  tt/6 k)  and  the 
equilibrium  kinetic  energy  of  order  E^O  =  ?r/2)  for  a 
wide  range  of  angles  0  . 


and 


The  mean  kinetic  energy  (4.42)  increases  sharply  with 
decreasing  0,  and  laser  cooling  is  possible  onlv  for  0  > 

(2  /I)'1'. 

One  can  see  that  the  minimal  equilibrium  energy  is 
achieved  for  a  small  angle  0  satisfying  the  condition 
(4.40).  From  Eq.  (4.33)  one  has 


is  the  mean  equilibrium  energy  reached  for  the  laser 
beams  with  orthogonal  polarizations.  Taking  into  ac¬ 
count  Eq.  (4. 38)  one  can  see  that  the  energy  E'™"  ob¬ 
tainable  for  small  angles  can  be  significantly  smaller  than 
£«,(#=  $)■ 

The  results  presented  above  for  the  equilibrium  en¬ 
ergy  distribution  have  been  derived  using  the  Fokker- 
Planck  equation  (3.18)  with  the  drift  and  diffusion  coef¬ 
ficients  averaged  over  the  optical  wavelength.  However, 
one  can  deduce  from  Eq.  (4.44)  that,  for  small  0,  the 
energy  E'”1"  is  smaller  than  the  maximum  of  the  effec¬ 
tive  potential  Uq .  Consequently,  an  atom  having  this 
energy  must  be  strongly  localized  [17],  This  means  that, 
rigorously  speaking,  the  equilibrium  energy  distribution 
l !'(:,]>)  should  strongly  depend  on  the  center-of-mass  po¬ 
sition  c,  and  that  localization  effects  must  play  a  signifi¬ 
cant  role  in  the  calculation  of  E«,.  Although  the  results 
obtained  using  the  averaged  Fokker- Planck  equation  are 
quite  similar  for  Je  =  1/2  and  3/2  atoms,  the  localization 
effects  make  these  two  systems  very  different.  It  has  been 
shown  earlier  that  atomic  motion  in  the  vicinities  of  the 
laser  field  nodes  is  responsible  for  a  dramatic  increase  of 
the  frictional  force  and  decrease  of  the  sub-Doppler-limit. 
monn  ntum  diffusion  that  may  eventually  lead  to  mean 


V.  CONCLUSIONS 

We  have  seen  that  in  the  one-dimensional  field  con¬ 
figuration  with  almost  parallel  linear  polarizations  of  the 
counterpropagating  laser  beams  the  considerable  increase 
of  the  effective  optical  relaxation  time  and  the  popula¬ 
tion  difference  gradient  near  the  nodes  of  the  field  may 
lead  to  significant  cooling  and  localization  of  the  atomic 
particles  having  s  Identical  effects  would  occur 
if  the  linearly  polarized  fields  were  replaced  by  two  cir¬ 
cularly  polarized  standing-wave  fields  (polarizations  e+ 
and  e_,  relative  phase  shift  0  <s(  1).  Although  the  cal¬ 
culation  has  been  limited  to  the  J9  —  -J  ground  state, 
one  would  expect  a  similar  qualitative  behavior  for  other 
values  of  Jg.  To  account  for  the  rapid  spatial  variation 
of  the  ground-state  population  difference,  new  analytical 
approaches  are  needed  since  methods  based  on  Fourier 
series  or  expansions  about  v  =  0  may  converge  very 
slowly.  The  need  for  new  analytical  techniques  is  under¬ 
scored  if  the  results  are  generalized  to  two-dimensional 
cooling.  In  this  case,  for  certain  field  configurations,  we 
find  a  dependence  of  the  friction  force  which  varies  as 
E j-  y  ~  nrr,y  !n(o2  +  o2),  where  ar  y  =  kvx  yrp. 

As  an  example  of  a  method  by  which  one  can  measure 
directly  the  dependence  of  the  spatially  averaged  cool¬ 
ing  force  on  angle  0,  we  should  like  to  cite  the  recent 
experiment  of  Grynberg,  Vallet,  and  Pinard  [18].  They 
measured  the  changes  in  field  intensity  as  two  copropa¬ 
gating  waves  traverse  a  medium  of  J  =  \  ground-state 
atoms  as  a  function  of  the  frequency  difference  6'  be¬ 
tween  the  waves.  The  changes  in  field  intensity  can  be 
related  dire  ’  ly  to  the  spatially  averaged  friction  force 
that  would  appear  in  the  sub-Doppler-limit  laser  cooling. 
The  two  fields  they  used  were  linearly  polarized  with  an 
angle  0  =  zr/ 2.  The  characteristic  width  they  found  in 
their  experiment  is  consistent  with  the  capture  range  that 
would  appear  in  the  sub-Doppler-limit  laser  cooling.  In 
the  similar  type  of  meiisurement,  by  changing  angle  0  be¬ 
tween  the  field  polarization  directions  one  might  see  the 
increase  of  the  signal  gradient  abo  t  S'  =  0  correspond¬ 
ing  to  the  increase  of  the  cooling  lorce  for  small  angles 
found  in  this  paper. 
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APPENDIX:  DERIVATION  OF  EQS.  (3.2) 
AND  (3.3) 

In  tliis  paper  we  consider  atoms  having  ground  and 
excited-state  angular  momenta  Jg  =  ^  and  Je  =  ^  or  ^ . 
The  atomic  density  matrix  sat  isfies  the  equation 

Jl  =  l  [/'•//]  +Kel(/»).  (AI) 

where  the  Hamiltonian  //  is  given  hy 


II  =  +  II M  +  //af. 


with  p  =  M  v  and  M  denoting  the  atomic  moment  uni  and 
mass,  respectively.  In  the  rotating-wave  approximation, 
the  atomic  Hamiltonian  //al  is  given  by 

Hai  =  -hb  ^  \Jc,iue)(Jt,wt\.  (A3) 

For  Je  =  the  interaction  part  of  the  Hamiltonian  H* r. 
which  describes  atom-field  coupling,  has  the  form 


-  hfy  g  COS  (2  +  *0  ^1/2)^-|/-’l 


COS  \2  +  J  ^ -t/'-OO/i/'-'l  •  (A') 


where  r  is  the  atomic  center-of-mass  position  In  the  case 
J,  =  t;  one  has 


-  -J  cos  ^1'  :t/L'X.'/i/-’l  +  -j\<  t/'-’l^  -  c<»s  ^~7j  +  k: 


I'  l/’.'l  +  -^=||  _  1  /■_>)(;/ 1  / 


'Fhe  relaxation  due  to  spontaneous-emission  processes  is  included  in  the  term  Rel(/<). 

I  sing  the  weak-field  condition  (1.1)  one  can  adial.at  a  ally  eliminate  the  excited-state  populations  and  excited 
ground-state  coherences  [1].  for  the  density  matrix  in  the  Wigner  representation  defined  In 


p(z.p)  =  (2  t/i)  1  J  p<;>  + 

in  the  =  7;  case,  one  finds  the  quasist  at  ionary  solution 


F'±i/-'  ~  ^  ^  /  F/2) vT1/-^'  ~  hhk.  ;)<  {±  jH  +  k  1 


+  /'.v?  1  /-’  .v  1  / -’ ( /'  +  :)' 


-il  T  Iv  +  k:  11 


1  /- 

P'±i/-"±i /-'(/'• -)  ~  +  l'-VvTi/:-yTi/L'(/'.-)cos(±0  -+-  2Lc) 


+/'vti/-’  st !/-’(/'  +  l>k.  :)  +  /'VT i/-_.  vt 1 /->(/'  -  *)]  +  1 


The  ground-state  populations  evolve  a 


(ll  P.vi  I/'-’  1  /'-'(/'■  -  )  -  ±  yjrl/'’  T  1/-  !/:’(/'  |/»F.  c  )<  "  *  -  M+  4  1  +  /I,  ^  1  /•_.  1  /•_.(/»  +  -j/iL,  c )(  11  *  "■**-']  ±  C.C 

1  f1 

+  2*  /  I^1  +,/'^vTi /'.'-T 1  /•-■(/'+  lik</.:)  +  (!  -  <lJ)l',±i/j,n/->(p+ hk<i.  :)}(!, 1,  (AH) 


where  the  integral  term  accounts  lor  I  lie  n  populal ion  ol  those  levels  resulting  from  the  emission  of  circularly  and 

linearly  polarized  photons  . .  the  excited  stale  [|  |]  Subsl  it  ul  mg  Fqs.  ( A  7 )  an. I  (  A  ,s  )  into  Fq  (AH)  for  the  ground- 

state  sublevel  populations.  Olie  arrives  at  the  (dosed  Set  of  equal  iolls 
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Accesion  For 


±  i  Ul 

,itp±l>'')  6(<s-'  +  r-'/i) 


2p±(P>  ~ )  +  cos(  ^ 0  +  2k: ) [/>±  ( /r  +  hk ,  :)  +  (>±{p  —  Ilk ,  : )] 


2  b 


+—  siu(T^  +  2k:)[f>±(ii  +  fik,:)  -  p±(p  -  hk,  .)] 


NTIS  CRA&I 
DTIG  TAB 

UnannOLnceo 
JuStif  IC3t>C 


+  - 


r/- 


By 


J  {( 1  +  >r)[p± (p  +  hk(q  +  1).  :)  +  p±(p  +  hk(<j  -  1).  :) 


+  l 

+2(>±(]>  +  hkii,:)cos(T0  4-  2k:)}  +  ( 1  -  <T)[/'t(p  + 1 2 * * * * 7 * * 10  11 12  +  1  )•  -)  , 

+Pt(p  +  hk((i  -  \),z)  +  2pT(p+  hk<i,:)cos(±0  A-2k: )]}</</, 

w hero  / >±  =  pg±i/-j  g±ipj- 

For  an  atom  having  ■],  —  4  t  he  corresponding  equations  are  v 


butio'  / 


(A  10) 


~(>±(lh  =  )  =  ~2^~frvT)  +  l;l  cos( ±0  +  +  rosfrO  +  ‘2A-r  )][/)*(/<  +  hk,:)  +  p±(p  -  hk.  c )] 

+  —  [2sin(±0  +  2k:)  +  si u(=pfl  +  2A-;)][/>-t(/>  +  hk.  ;)  -  />±(/<  -  hk.  i)]^ 

+  J  1  +  '/-){V±(/>+  hk(<i  +  1).  :)  +  .r)/)±(;,  +  hk(<]  -  1),  c) 

+P±ii>  +  hki/,  j)[‘Jcos(±d  +  2Ac)  +  cos(^P  I  2A;)]} 

+2(  \  -  <r)[t>T(p+  hk{<i  +  H.c)  +  /'T(/»+  hk(<i  -  1).:) 

+2/>t(/>  +  /)<■</.  : )  cos (iff  +  2k:  )]}(/</.  (All) 

Equations  (A  10)  and  (All)  are  exact  under  the  uvak-lield  assumption  (I  I)  If,  in  addition,  one  assumes  that  an 
atomic  momentum  is  much  larger  than  that  of  the  photon,  that  is 


P  »  hk. 

one  can  expand  these  equations  to  second  order  in  hk  to  arrive  at  Kqs.  (2.2)  ami  (2.2)  for 


(A  12) 

and  4-  respectively. 
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